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Abstract 

Motivated by the recent works of Refs. [U |2] where a model of inflation has been 
suggested with non- minimally coupled massive vector fields, we generalize their work 
to the study of the bouncing solution. So we consider a massive vector field, which 
is non-minimally coupled to gravity. Also we consider non-minimal coupling of vector 
field to the scalar curvature. Then we reconstruct this model in the light of three forms 
of parametrization for dynamical dark energy. Finally we simply plot reconstructed 
physical quantities in flat universe. 
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1 Introduction 

Nowadays it is strongly believed that the universe is experiencing an accelerated expansion. 
The observation data confirm it such as type la supernovae [3] in associated with large scale 
structure [1] and Cosmic Microwave Background anisotropics [5] have provided main evi- 
dence for this cosmic acceleration. In order to explain why the cosmic acceleration happens, 
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many theories have been proposed. The standard cosmological model (SCM) furnishes an 
accurate description of the evolution of the universe, in spite of its success, the SCM suf- 
fers from a series of problems such as the initial singularity, the cosmological horizon, the 
flatness problem, the baryon asymmetry and the nature of dark energy and dark matter, 
although inflation partially or totally answers some of these problems. Inflation theory was 
first proposed by Guth in 1981 [B]. Inflation is a period of accelerated expansion in the early 
universe, it occurs when the energy density of the universe is dominated by the potential 
energy of some scalar field called inflaton. Currently all successful inflationary scenarios are 
based on the use of weakly interaction scalar fields. Scalar fields naturally arise in particle 
physics including string theory and these can act as candidates for dark energy. So far a wide 
verity of scalar field dark energy models have been proposed. These include quintessence [7], 
K-essence [S], tachyon [H], phantoms [TU], ghost condensates [TT] and so forth. Two main 
reason for use of scalar fields to explain inflation are natural homogeneity and isotropy of 
such fields and its ability to imitate a slowly decaying cosmological constant [1]. However, no 
scalar field has ever been observed, and designing models by using unobserved scalar fields 
undermine their predictability and falsifiability, despite the recent precision data. The latest 
theoretical developments (string landscape) offer too much freedom for model-building, so 
higher spin fields generically induce a spatial anisotropy and the effective mass of such fields 
usually of the order of the Hubble scale and the slow-roll inflation dose not occurs [12] . 
Then an immediate question is, can we do Cosmology without scalar fields? The authors 
of [H [2] have shown that a successful vector inflation can be simultaneously surmounted in 
a natural way, and isotropy of the vector field condensate be achieved either in the case of 
triplet of mutually orthogonal vector field [13] . In spite of inflation success in explaining the 
present state of the universe, it dose not solve the crucial problem of the initial singularity 
[13] . The existence of an initial singularity is disturbing, because the space-time description 
breaks down "there". Non-singular universes have been recurrently present in the scientific 
literature. Bouncing model is one of them that was first proposed by Novello and Salim 
[T^ and Mlnikov and Orlov [TB] in the late 70's. At the end of the 90's the discovery of 
the acceleration of the universe brought back to the front the idea that p + "ip could be 
negative, which is precisely one of the conditions needed for cosmological bounce in OR, 
and contributed to the revival of nonsingular universes. Bouncing universe are those that 
go from an era of acceleration collapse to an expanding era without displaying a singularity 
[T7] . Necessary conditions required for a successful bounce during the contracting phase, the 
scale factor a{t) is decreasing , i.e. d < 0, and in the expanding phase we have d > 0. At 
the bouncing point, d = 0, and around this point d > for a period of time. Equivalently 
in the bouncing cosmology the Hubble parameter H runs across zero from < to > 
and = at the bouncing point. A successful bounce requires around this point. 
The remainder of the paper is as follows. In section 2 and 3, we will consider vector field 
action where proposed in Refs. [H [2] and study bouncing solution of this model. In section 
4 and 5, we will reconstruct physical quantities for this model and also will plot the corre- 
sponding graphs. Finally we will apply three parametrization and compare them for this 
model. 
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2 Vector field foundation 

We consider a massive vector field, which is non-minimally coupled to gravity, [H [2] . The 
action is given by 

S = J d'x^g {j^R - \f,.F^'' - Im'U.U^ + l^RU.U''^ , (2.1) 

where F^^y = d^Uy — dpU^^, and ^ is a dimensionless parameter for non-minimal coupling. 
We note that, the non-minimal coupling of vector field is same with conformal coupling of a 
scalar field in case ^ = 1/6. We adopt FRW universe with the metric signature of ( — h ++). 
The equations of motion are given by 



RiJiv — -^R9tiu — 8vrG 



(2.2) 



VyF""^ - m^W + ^RW = 0. (2.3) 

Where the right hand side of equation (12.21) is the energy- momentum tensor of the vector 
field Ui. The variation of the action with respect to Ui yields the following equations of 
motion, 

^V't/o - ^d^ili - m^Uo + iRUo = 0, (2.4) 

Ui + -{Ui - d^Uo) - diUi + hdiidkUk) - V^Ui) + m'^Ui - iRU = 0, (2.5) 
a 

Where a is the scale factor, the dot denotes the derivative with respect to the cosmic time 
and the summation over repeated spatial indices is satisfied. By considering the quasi- 
homogeneous vector field {diUa = 0) and Eq. (12. 4p imply Uq = 0, so that from Eq. (12. 5 p we 
obtain 

tli + HUi - Q^{H + 2H^ + ^)Ui + m^Ui = 0. (2.6) 
By using acceleration relation ^ = — ^y^(p + Sp) we achieve as, 

H + H' = ^ (2Uf - 4(1 + QmuA + m^H' - m'Ul - -^iui) , (2.7) 

where H = ^, R = 6(if ^ + 2if ^ + ^) and Rq = H + H"^ are Hubble 's parameter, Ricci scalar 
and first component of Ricci tensor, respectively. As we know, a dynamical vector field has 
generally a preferred direction, and to introduce such a vector field may not be consistent 
with the isotropy of the universe. In fact, the energy- momentum tensor of the vector field f/^ 
has anisotropic components. However, the anisotropic part of the energy-momentum tensor 
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can be eliminated by introducing a triplet of mutually orthogonal vector fields. In that case, 
we obtain the energy density p and the pressure p of the vector fields 



P = — 



lui" - 3(1 + QOHUiU, + g^f/'i/' + ^m^Uf 



,2 ^« 



(2.8) 



1 



3 3A;f 
- 3(1 + Qi)HU,U, + %U^H^ - -m'^Uf + -^Uf 



1-^ 



(2.9) 



Now to introduce a change of variable 0j = — (for more detail see Ref. [1]), equation (12. 4p 
make change to, 



i + m'<\) + + (1 - ^i){H + 



6^k 



0. 



(2.10) 



Then we consider ^ = 1/6 and obtain the basic equations of motion for a curved universe in 
terms of 0i, 

4>i + 3H^+(m^--^]4>i = 0, (2.11) 



(2.12) 



H + H^ = -47rG(202 _ m''(Pl), (2.13) 

One can see where equations of motion of vector field is reduced to minimally coupled massive 
scalar fields. So energy density p and the pressure p for the vector fields are derived in terms 
of (pi in case = 1/6 in the form. 



3;2 3 2i2 ^^i2 



-m 
2 



2a? 



3 'r, 3 9,2 ^ l1 



(2.14) 
(2.15) 



Now we are going to consider behavior of the different values of parameter ^ for vector field. 
We solve numerically Eq. (12. 6p for K = 0,+l,— 1 which implies the fiat, close and open 
universe respectively. The Fig.l shows graph of the vector field with respect to time in all of 
cases K. One can see where vector field has oscillation behavior and the magnitude slowly 
decrease with respect to time evolution. Also by increasing the parameter ^, the magnitude 
of vector field will increase, but the period of oscillation is constant. We note that negative 
values of ^ actually is the same of above result. 

As above mention we suggest following solution for Ui{t) 

Ui (t) = VAe~^^ cos{mt + 9), (2.16) 

where the parameter A describes the oscillating amplitude of the field with dimension of 
[mass]^. Also A is relation with the parameter i^, this solution implies the damping magni- 
tude of the oscillating vector field. 
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Figure 1: Graphs of vector fields in term of time. The sohd, dash and doted hnes represent 
^ = 1 ,| and 0.5 respectively. 

3 Bouncing behavior 

We will start with a detailed examination on the necessary conditions required for a successful 
bounce. During the contracting phase, the scale factor a{t) is decreasing, i.e., d < 0, and in 
the expanding phase we have d > 0. At the bouncing point d = and around this point 
d > for a period of time. Equivalently in the bouncing cosmology the Hubble parameter 
H runs across zero from H < to H > and if = at the bouncing point. A successful 
bounce requires around this point 

H = -A7tG{p + p) + -^>0. (3.1) 
At the point where the bounce occurs, Eqs. ( 12.81) and (12. 9p reduce to 

P. = ^(t>.' + m^[/f)-^et/f, (3.2) 

P.= ^(f>/-^^t/^) + ^ef/^ (3.3) 
On the other hand, a successful bounce from Eqs. (12.60 . (12. 7p and ( 13. ip obtain in the form, 

U- < \m'Ul + ^^Uf. (3.4) 

This result is similar to slow roll inflation. This means that one requires a fiat potential 
where give rise to a point bounce for the model of vector field. From conditions (13. ip . (13.40 
it is clear that if we have bouncing solutions in open universe, then we have such behaviour 
for fiat and closed universe as well. Now we solve above equation numerically by different 
value of ^ on the curved universe that is plotted in Fig. 2. 
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Figure 2: The graphs of the Hubble's parameter for ^ = 1/6, AnG = 1, m = 1 and k = 
0, +1, -1 by choosing 0(0) = 1, 0(0) = 0.1, a(0) = 1 and H{0) = 0.01. 
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Figure 3: The graphs of the scale factor for ,^ = 1/6, AivG = 1, m = 1 and k = 0, +1, —1 by 
choosing 0(0) = 1, 0(0) = 0.1, a(0) = 1 and H{0) = 0.01. 

One can see the Hubble parameter H running across zero in any three cases of k. In all 
cases of k, we have H < to H > where implies to go from collapse era to an expanding 
era, and this result will not change for the different values ^ in all of k. Also in Fig.3, we 
can see the behaviour of scale factor in terms of time for different values of k. It is clear 
that during the contracting phase, the scale factor a{t) is decreasing, i.e., a < 0, and in the 
expanding phase we have a > 0, so the point where a = is bouncing point. 
Therefore, in the vector field dominated universe we have a successful bouncing point in close 
and flat universe but a turn-around point in open universe. The bounce can be attributed 
to the negative-energy matter, which dominates at small values of a and create a significant 
enough repulsive force so that a big crunch is avoided. 
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4 Reconstruction 



Now we are going to present a reconstruction process for vector field in the curved universe 
by ^ = 1/6 . In this section, potential and kinetic energy are reconstructed with respect to 
redshirt z. Also we obtain the EoS in term of z. After that three type parametrization are 
represented for the EoS. By using it we consider cosmology solutions such as the Eos, the 
deceleration parameter and vector field. The stability condition of this system is described 
by quantity of the sound speed. We rewrite Eqs. (12.141) and (I2.15P in term of the effective 
potential energy V and the effective kinetic energy K as the following form. 



3 3A; ^ ^ 

P=2* + - = S'f + 3f . (4-1) 



p.|^.^-|,„V?H-^^?^3i--3V'-±, (4.2, 
p + p = Qk--^. (4.3) 



Then we can write the Friedmann equations as following 

3M^{H^ + \) = p,^ + p = p^ + 3k + 3V, (4.4) 

2Ml{H -^) = -pra- P-P=-Pra-Qk+^, (4.5) 

where p^ is the energy density of dust matter. Also from Eqs. (14.11) and (14. 2p . we obtain 
relationship between the Eos with V and K as, 

p 3k-3V-^ 2--!^ 

uj = - = . ^ = -1 + 4.6 

P 3K + 3V 1 + ^ 

We obviously have 



V + 3k> 


k 

3^ ~ 


^ u > 


-1, 


V + 3k < 


k 

3^ ~ 


^ id < 


-1, 


v + k = 


k 

3^ ~ 


=^ u = 


-1, 



(4.7) 



By using Eqs. (14.41) and (14.51) we can write 



K^Z£Il.^^H-±) + -lL, (4, 



V = ^i^H' + H+^)-^--^. (4.9) 
3 6 6a^ 
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As in the present model, the dark energy fluid does not couple to the background fluid, the 
expression of the energy density of dust matter in respect of redshift z is 



p^ = 'iMlEl'a^,{\^z)\ (4.10) 

where ^raa is the ratio density parameter of matter fluid and the subscript indicates the 
present value of the corresponding quantity. By using the equation Wz = ^ (ao is quantity 
given at the present epoch) and its differential form in following have, 



To introduce a new variable r as. 



(4.12) 



we rewrite the equation of motion of vector field against z as, 

2r(l + zfV'l + 2r(l + ^)(1 + Hl)U[ - r'(l + zfV[ 



2m\, 2k 



+r'(l + z)Ui - rUi + —^U, - ^-^ (1 + zfUi = 0, (4.13) 

K, V can be rewrite as following 

K = -\MlHlnU^ + zf + \MlHlr\l + z) + ^(1 + z)\ (4.14) 

2k 1 1 k 

V = MlHlr + —{l + zf- -MlHlil + zY - -MlHln^o{l + " + ^f- (4-15) 

By using Eqs. (14. 6p . (14. 8 p and (14. 9 p we obtain following expression for the EoS, 

(l + ^V-3r + ^(M,^-2) 
"3r-3a^o(l + ^)^ + ^(M^2 + 2)- 

Then we obtain following equation for r{z) 

r{z)=n^^{l + zf + {l-Q^o)e''^'^ 
+ao f [uj{~z){2 + Ml) + (2 - Ml)] (1 + ~z)e-''^'^ d~z. (4.17) 



where (5{z) = /; ^dz and ao = ^x^- 



r{z) = fi^o(l + z)' + (1 - fi„o)e3/o i±^ci5. (4.18) 
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Also we have following expression for deceleration parameter q 



H (l + z)r'-2r 
^(-) = -1 - ^ = ^ h • (4.19) 



Now we consider the stability of this model by use the hydrodynamic analogy and judge on 
stability by examining the value of the sound speed. Of course this is a simple approach, the 
perturbations in vector inflation are much richer than in hydrodynamic model, see recent 
interesting works in [121 [20]. The sound speed can be obtained by the following equation, 

' P' -9fi„o(l + ^)2 + 3r' + 2^|^(M2 + 2)' 

in order to deal the stability of our model, the sound speed must become > 0, so we can 
obtain from above equation following condition 

r(z)>c.„o(l + .)^-^^^. (4.21) 



5 Parametrization 

Now we consider the three different forms of parametrization as following and compare them 
together. 

Parametrization 1: First Parametrization has proposed by Chevallier and Polarski [21] .and 
Linder [22], where the EoS of dark energy in term of redshift z is given by, 

u{z)=uJo + ^^^. (5.1) 
1 + z 

Parametrization 2: Another the EoS in term of redshift z has proposed by Jassal, Bagla 
and Padmanabhan [22] as, 

u{z) = uo + (5.2) 
(1 + z)^ 

Parametrization 3: Third parametrization has proposed by Alam, Sahni and Starobinsky 
[23]. They take expression of r in term of z as foUowoing, 

r{z) = fi^o(l + + Ao + Ai(l + z) + A2{1 + zf. (5.3) 

By using the results of Refs. [25], [26], [2Z] and [21], we get coefficients of parametrization 1 
as fimo = 0.29, Uq = —1.07 and Ua = 0.85, coefficients of parametrization 2 as Qmo = 0.28, 
ujq = —1.37 and Ub = 3.39 and coefficients of parametrization 3 as flmo = 0.30, Aq = 1, 
Al = —0.48 and A2 = 0.25. The evolution of u{z) and q{z) are plotted in Fig. 3. Also, 
using Eqs. f l4.14p and fl4.15p and the three parametrization, the evolutions of K{z) and V{z) 
are shown in Fig. 5 and Fig. 6 respectively. We note that graphs simply represent only in the 
fiat universe. From Figs. 4, 5 and 6, we can see parametrization 1 and 3 are same nearly and 
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Figure 4: Graphs of the EoS and deceleration parameter in respect of redshift z. The sohd, 
dot and dash hnes represent parametrization 1, 2 and 3, respectively. 



3-" 




Figure 5: Graphs of the reconstructed K in respect of redshift z. The solid, dot and dash 
lines represent parametrization 1, 2 and 3, respectively. 

have slightly different from parametrization 2. The EoS for any parametrization show in Fig. 
(4) so that they running cross to —1. Acceleration for all of parametrization shows to tend to 
the positive value. The K and V increase for parametrization 1 and 3, but in parametriza- 
tion 2 increase (decrease) for the V [K\ One can see that parametrization 1 and 3 satisfy 
condition + > and parametrization 2 satisfy condition K = Constant. In order 
to by Eq. (14. 7p . we have a; > — 1 (a; = — 1) when parametrization 1 and 3 (parametrization 
2). This is mean that parametrization 2 is better than others parametrization. 

In Fig. 6, we can see the variation of the vector field against redshift z. One is obviously 
showed slightly difference between all of parametrization. 
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Figure 6: Graphs of the reconstructed V in respect of redshift z. The sohd, dot and dash 
hues represent parametrization 1, 2 and 3, respectively. 

6 Conclusion 

In this paper, we have studied the bouncing solution in curved universe which proposed by 
the model of a massive vector field, Ui, non- minimally coupled to gravity. For our purpose 
we have derived the corresponding energy density, pressure and Friedmann equation for this 
model. Also we have obtained the bouncing condition as Eq. fl3.4p . From this condition, and 
also essential condition (13. ip . it is clear that if we have bouncing solutions in open universe, 
then we have such behaviour for fiat and closed universe as well. After we plot the Hubble 
parameter in term of time by figures 2, for different k, we understood that our model predict 
the bouncing behavior for all cases of k. Fro these figures one can see that the Hubble 
parameter H running across zero in any three cases of A;. In all cases oi k, we have H < 
to if > where implies to go from collapse era to an expanding era, and this result will not 
change for the different values ^ in all of k. After that in figures 3 we have shown that during 
the contracting phase, the scale factor a{t) is decreasing, i.e., a < 0, and in the expanding 
phase we have a > 0, so the point where a = is bouncing point, and this figure is consistent 
with the results of Fig 2. 

After that we have investigated an interesting method as the reconstruction of the non- 
minimally coupled massive vector field model with the action (12. ip . Our aim was to see 
whether the non-minimal coupling vector field can actually reproduce required values of ob- 
servable cosmology, such as evolution of the EoS and the deceleration parameter in respect to 
the redshift z. We have reconstructed our model in three different forms of parametrization 
for massive vector field. In Fig. 4 we have found the EoS crossing —1 in all of parametriza- 
tion. The variation of reconstructed kinetic and potential energy against z have plotted in 
Figs. 5 and 6, where the parametrization 2 in addition is better than two others parametriza- 
tion because K + V = Constant. Also we have investigated the stability of this system and 
have obtained a condition by the sound speed in all of curvatures. Finally we note that 
reconstructed physical quantities have just executed in fiat universe and one is suggested for 
open and close universe as future work. 
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z 

Figure 7: Graphs of the reconstructed f/j in respect of redshift z. The sohd, dot and dash 
hues represent parametrization 1, 2 and 3, respectively. 



7 Acknowledgment 

The authors are indebted to the anonymous referee for his/her comments that improved the 
paper drastically. 



References 

[1] A. Golovnev, V. Mukhanov and V. Vanchurin, JCAP 0806, 009 (2008). 
[2] T. Chiba JarXiv:0805.4660i ^3 [gr-qc]2008 

[3] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) 
astro-ph / 9805201J ; S. Perlmutter et al. [Supernova Cosmology Project Collaboration], 



Astrophys. J. 517, 565 (1999) |astro-ph/9812133| ; P. Astier et al, Astron. Astrophys. 
447, 31 (2006) |astro-ph/0510447[ ; 



[4] K. Abazajian et al. [SDSS Collaboration], Astron. J. 128, 502 (2004) astro-ph/0403325 



K. Abazajian et al. [SDSS Collaboration], Astron. J. 129, 1755 (2005) |ast ro-ph/04102 39|. 
[5] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 



|a stro-ph/0302209] ; D. N. Spergel et al., |astro-ph/0603449, 
[6] A. H. Guth, Phys. Rev. D 23, 347, (1981). 



12 



[7] Peebles, P.J.E., Ratra, B.: Astrophys. J. 325, L17 (1988); Ratra, B., Peebles, P.J.E.: 
Phys. Rev. D 37, 3406 (1988); Wetterich, C: Nucl. Phys. B 302, 66 8 (1988); C aldwell, 
R.R., Dave, R., Steinhardt, P.J.: Phys. Rev. Lett. 80, 1582 (1998). 'astro-ph/9708069; 
Liddle, A.R., Scherrer, R.J.: Phys. Rev. D 59, 023509 (1999). astro-ph/9809272i Zlatev, 
I., Wang, L.M., Steinhardt, P.J.: Phys. Rev. Lett. 82, 896 (1999). astro-ph/9807002; 
Faraoni, V. Phys. Rev. D62, 023504 (2000); Huang, Z.G., Lu, H.Q., Fang, W.: Class. 



Quantum Gravity 23, 6215 (2006). |ihep-th/0604160 



[8] Armendariz-Picon, C, Mukhanov, V.F., Steinhardt, P.J.: Phys. Rev. Lett. 85, 4438 
(2000). |astro-ph/ 0004134[ Armendariz-Picon, C, Mukhanov, V.F., Steinhardt, P.J.: 
Phys. Rev. D 63, 103510 (2001). astro -ph/ 0006373} Chiba, T., Okabe, T., Yamaguchi, 
M.: Phys. Rev. D 62, 023511 (2000) 

[9] Sen, A.: J. High Energy Phys. 0207, 065 (2002); Padmanabhan, T. Phys. Rev. D 66, 
021301 (2002).; Setare, M. R. Phys. Lett. B 653, 116 (2007); M. R. Setare, J. Sadeghi, 
A. R. Amani, Phys. Lett. B 673, 241, (2009). 



[10] Caldwell, R. R. Phys. Lett. B 545, 23 (2002). |astro-p h/9908168; Caldwell, R. 
R., Kamionkowski, M., Weinberg, N. N. Phys. Rev. Lett. 91, 071301 (2003). 



astro-ph/ 0302506, Nojiri, S., Odintsov, S.D.: Phys. Lett. B 562, 147 (2003). 



hep-th/0303117; Nojiri, S., Odintsov, S. D. Phys. Lett. B 565, 1 (2003). hep-th/0304131 



Faraoni, V. Class. Quantum Gravity 22, 3235 (2005); Setare, M. R. Eur. Phys. J. C 50, 
991 (2007) 



[11] Piazza, F., Tsujikawa, S., J. Cosmol. Astropart. Phys. 0407, 004 (2004). |hep^th/0405054 
[12] L. H. Ford, Phys. Rev. D 40, 967 (1989). 
[13] Armendariz-Picon C, JCAP 07(2004)007. 



[14] A. Borde, A. Guth, A. Vilenkin, Phys. Rev. Lett. 90, 151301 (2003), |gr^/0110012 
[15] M. Novello and J. Salim, Phys. Rev. D 20, 377 (1979). 
[16] V. N. Melnikov and S. V. Orlov, Phys. Lett. A 70, 263 (1979). 
[17] M. Novello, S. E. Perez Berdiaffa. E-Xiv:0802. 16341 [astro-ph]. 

[18] S. Zhang, B. Chen, Phys. Lett. B 669, 4 (2008); M. R. Setare, J. Sadeghi, A. Banijamah, 
Phys. Lett. B 669, 9, (2008); M. R. Setare , J. Sadeghi and A. R. Amani ,Int. J. Mod. 
Phys. D 18, 1291, (2009). 

[19] B. Himmetoglu, C. R. Contaldi, M. Peloso, Phys. Rev. D 79, 063517, (2009); B. Him- 
metoglu, C. R. Contaldi, M. Peloso. [arXiv:0909.3524 yl [astro-ph. CO]; B. Himmetoglu, 
C. R. Contaldi, M. Peloso, Phys. Rev. Lett. 102, 111301, (2009). 

[20] A. Golovnev. rarXiv:0910.0173v 4 [astro-ph. CO]; A. Golovnev, V. Vanchurin, Phys. Rev. 
D 79, 103524, (2009). 



13 



[21] M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 213 (2001). 
[22] E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003). 

[23] H. K. Jassal, J. S. Bagla and T. Padmanabhan, Mon. Not. Roy. Astron. Soc. 356, Lll 
(2005). 11 

[24] U. Alam, V. Sahni, T. D. Saini and A. A. Starobinsky, Mon. Not. Roy. Astron. Soc. 
354, 275 (2004). 

[25] A. G. Riess et al. [Supernova Search Team Collaboration], [arXiv:astro-ph/0611572 . 

[26] Y. G. Gong and A. Z. Wang, Phys. Rev. D 75, 043520, (2007). 

[27] Y. Wang and P. Mukherjee, Astrophys. J. 650, 1 (2006). 

[28] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005). 



14 



